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The non-Markovianity is a prominent concept of the dynamics of the open quantum systems, 
which is of fundamental importance in quantum mechanics and quantum information. Despite of 
lots of efforts, the experimentally measuring of non-Markovianity of an open system is still limited to 
very small systems. Presently, it is still impossible to experimentally quantify the non-Markovianity 
of high dimension systems with the widely used Breuer-Laine-Piilo (BLP) trace distance measure. In 
this paper, we propose a method, combining experimental measurements and numerical calculations, 
that allow quantifying the non-Markovianity of a N dimension system only scaled as N 2 , successfully 
avoid the exponential scaling with the dimension of the open system in the current method. After 
the benchmark with a two-dimension open system, we demonstrate the method in quantifying the 
non-Markovanity of a high dimension open quantum random walk system. 

PACS numbers: 03.65.Yz, 42.50.-p, 03.67.-a 


I. INTRODUCTION 

No real physical systems could be regarded as a purely 
closed system, as they are inevitably interacting with 
environments. Therefore, the dynamic of the open 
quantum systems is in the central of the fundamental 
quantum mechanics and quantum information science, [lj 
Non-Markovianity [2} is the prominent concept in open 
systems and attracts a lot of attention from theoretical 
and experimental aspects. M3 It has been shown that 
the non-Markovianity can be exploited as useful resource 
in quantum technology. For examples: it may be used 
to improve efficiency of quantum information processing 
and communication; [4] it is benefit in quantum metrol¬ 
ogy; i and it may be used to improve the security in 
continuous-variable quantum key distribution, 0 etc. It 
is therefore important to quantify the non-Markovianity 
of a quantum system. However, despite of the extensive 
investigation, characterization and measurement of the 
non-Markovianity of an open system is still limited to 
very small systems, I la llJl mostly one-qubit system in 
experiments, flij ] 

The measurement base on trace distance, proposed by 
Breuer, Laine and Piilo (BLP), [g] is one of the most 
popular definitions for non-Markovianity, and is widely 
used in theoretical and experiment investigations. |15l.ll6l| 
To quantify the non-Markovianity by BLP measurement, 
one needs to find a pair of initial states to maximal a 
function based on the trace distance. Only when the 
interactions between the system and the environment 
are exactly known (the dimension of the whole system 
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is not too huge), the optimal state pair and the measure 
of the non-Markovianity can be found numerically. For 
very limited quantum open systems which can be exactly 
solved, such as the Jaynes-Cummings model 0 and 
quantum Brownian motion model, fl8j the measure of the 
non-Markovianity can be analytically found. However, 
for general quantum open systems where we have no 
exact information about the interaction between the 
system and the environment (or the dimension of the 
whole system is large), the optimal initial state pair can 
only be found experimentally by scanning the dynamics 
of the whole initial state spaces which is a tough task even 
for two-dimension systems. In a typical experiment [b9j 
to quantify the non-Markovianity in a two-dimension 
quantum open system, in order to achieve a reasonable 
accuracy, total 5000 states’ dynamics were measured. 
Even worse, the number of scanning states grows expo¬ 
nentially with the dimension of the system. For a system 
containing two spin-1/2 qubits (which can be viewed 
as a 4-dimension system), one need to scan about ZJ 4 
different initial states to obtain the non-Markovianity, 
where D is the number of samples for each degree of 
freedom. Typically, D should be 100 or even more to 
ensure the accuracy. Therefore even quantification the 
non-Markovianity of this simple two qubits system is 
beyond our current experimental ability. 

In this work, inspired by the idea of standard quantum 
process tomography, [13 we propose an efficient method, 
combining handful experiments (polynomial scaled with 
the dimension of the system) and the numerical optimiza¬ 
tion method to measure the non-Markovianity for high 
dimension systems (The systems containing more than 
one qubits can be regarded as high dimension systems). 
This method require no prior information about the in¬ 
teractions between the system and its environment. Due 
to the linearity of the dynamics of the quantum system, 
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we need only experimentally measure the dynamics of 
some linearly independent states of the system, and the 
dynamics of the whole state space can be rebuilt by linear 
combination of these experimental results. We then find 
the optimal state pair in quantifying non-Markovianity 
through numerical calculation based on these experiment 
data. The number of the measurement is scaled as N 2 , 
where N is the dimension of the quantum system. Using 
this method, the former intractable non-Markovianity 
measurement of high dimension system can be easily 
investigated. After the benchmark on a two-dimension 
open system, we demonstrate our algorithm to quantify 
the non-Markovianity of a high dimension open quantum 
random walk system. 


II. METHODS 

The BLP measurement of the non-Markovianity [fi] is 
defined on the trace distance of two state p±, P2, that 
is, D(p 1 ,p 2 ) = |tr|pi - p 2 \ where |A| = \fAAA. If A 
is Hermitian, tr|A| = JA |Aj| is the sum of the absolute 
value of all the eigenvalue of matrix A. This quantity 
describe the distinguishability between the states p\ and 
p 2 : if it is zero, the two states are indistinguishable, 
otherwise, they are distinguishable. Based on this 
definition, the measure of the non-Markovianity of an 
open system can be defined as, 

M = max pi 2(0) / dtS(t,pi,2{0)), (1) 

J( 5>0 

where <5(i, Pi, 2 (0)) = -^D[pi(t), p 2 (t)\ is the change rate 
of the trace distance. Pi{t), i= 1, 2 is the density 
matrix of the open system at time t with the initial 
state Pi(0). The time-integration is extended over all 
time intervals in which <5 is positive, and the maximum 
should be optimized over all pairs of initial states. 
Roughly speaking, the integral intervals stand for the 
time intervals when the information flows back to the 
system from the environment. 

The most difficult task to measure the non- 
Markovianity is to find a pair of states, pi, p 2 that 
maximize Eq. ©. Generally, it need experimentally 
scan the state pairs in the whole parameter space. 
Some simplifications can be made. [19j, [21[ It has been 
rigorously proven that the optimal states pair should 
be on the boundary of the physical state space and the 
states pair are orthogonal each other. Therefor, we can 
only scan the boundary of physical state space (For two- 
dimension case, these states are pure states). Another 
simplification to the measure was demonstrated in Ref. 
Il9l which illustrated that the measure can be obtained 
efficiently in an arbitrary neighborhood of any fixed state 
in the interior of the state space. That is, it needs only 
scan one state of the pair in the physical state space. 
This can dramatically reduce the experimental work. 
However, the number of the experiments is still too large 


and will exponentially increase with the dimension of 
the systems. Therefore, the non-Markovianity in higher 
dimension is still intractable to quantify with the current 
method. Now we introduce another scheme to simplify 
the experimental quantification of non-Markovianity in 
an open system which make the high dimension system 
reachable. 

Following the idea of the quantum process tomogra¬ 
phy, [20| the state of a quantum open system with N 
dimension can be expressed as NxN density matrix. 
Any density matrix can be expanded by N 2 linear 
independent bases, 

P x mn = (|m)(n| + |n)(m|)/2, (m > n) 

Prnn = K\ m )( n \ ~ l™)M)/ 2 , i™ > Tl) 

Pm. = NH)> (2) 

where \m) (to = 1,2,- •• , N) is the basis vector of 
the system. The operators p^ nn (p^n) play the similar 
role of the pauli matrices a x ( a v ) in two-dimension 
systems. Without loss of generality, we assume that 
the system and the environment is in a product state 
at the initial time t = 0, i.e. p( 0) = p s (0) <g> p e { 0), where 
p(t), p s {t ), p e (t) are the density matrices of the whole 
system (system+enviroment), the quantum system and 
the environment at time t, respectively. Using the above 
introduced bases, the state of the open system at any 
time t can be written as: 

Ps(t) = aX muP X mn{t) + OLW + a mPm(t) ,(3) 

m>n m>n m 

where a^ nn , a v mn , are time independent constants 
determined by the initial states. p l mn {t) i = x,y, 0 is 
the dynamics of the bases. It suggests that the dynamics 
of the open system can be completely determined by the 
dynamics of the bases. 

In experiments, the dynamics of the N 2 bases can be 
obtained by the following procedure: 

1. Prepare the initial states of the system to |to), 

where m=l, 2, ■■■, N. \m ) can be any set 

of complete and orthogonal vector bases of the 
system. By measuring the dynamics of the open 
system, we obtain Pm(t)- 

2. Prepare the initial states of the system to (| to) + 

\n))/y/2 ( to > n), measuring the dynamics of the 
open system. We obtain p x mn {t) + + P°n (*))• 

3. Prepare the initial states of the system to (|to) + 
i\n))/y/2 ( to > n), measuring the dynamics of the 
open system. We obtain — p^„(t) + |(p^(t) + 
P°n(t)). 

We therefore have the dynamics of the N 2 bases of the 
open system. Using the dynamics of these bases, the BLP 
measure of the non-Markovianity of the open system can 
be achieved by numerically optimizing the parameters, 
a mm a mn an d C, of the initial states by computer using 
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FIG. 1: (Color online) Comparing the trace distance directly 
measured from experiment and the results obtained from our 
method. The blue circles show the trace distance as function 
of time of the experimentally determined optimal state pair, 
by scanning the whole state space. The red squares show the 
trace distance of the optimal pair as function of time, obtained 
by four basis dynamics followed by numerical optimization, as 
described in the main text. 

Eq. © ,and Eq. Q through deepest descent algorithm. 
The nice scale of this method make it possible to apply 
for the high-dimension system which is intractable for 
the traditional method. It is worth noting that the 
basis introduced here is not unique. Any N 2 such linear 
independent states are enough for the procedure. The 
simplifications introduced in Ref. [U and Ref. |l9j, can 
also be used to reduce the numerical optimization efforts 
in this method. 


III. RESULTS AND DISCUSSION 
A. Benchmark for the two-dimension system 

To demonstrate the power of our scheme, we first 
make benchmark tests on the well studied two-dimension 
open system. In the typical experiment, [19] the two- 
dimension open quantum system is provided by the po¬ 
larization of a photon which coupled to the environment 
through its frequency degree. To quantify the non- 
Markovianity of this system, 5000 different states in the 
Bloch surface have been scanned to find the optimal 
pair in Eq.(l). The conclusion in j2l] has been used to 
simplify the experiment, which states that the optimal 
state pairs are orthogonal pure states, and therefore, 
function J\f in Eq.(l) only depends on the angles between 
two states on the Bloch surface. 

A 2 x 2 density matrix can be expanded by identity 
matrix I and pauli matrixes a x , a y , a z as p s = ^(I+a-a), 
where vector a is on the surface of the Bloch sphere. 
Using our method, non-Markovianity measure of this 
open system can be determined by the dynamics of the 
pauli matrices, i.e. where i=x, y, z which can 

obtained from the dynamics of the four initial states, 11), 



FIG. 2: (color online) Cartoon of the 3-step open quantum 
walk system, with X=2. The blue blocks denote the polarized 
beam splitter. This quantum walk system includes both the 
lattice (location of the photon) and the coin (polarization) 
which is a intrinsically high-dimension system. 


I- 1 ) > 73 ( 11 ) + I - !)) and 73GI)+*|- 1 )) (where | - 1) 
and |1) are two eigenvectors of the polarization photon). 

Fortunately, the dynamics of these four initial states 
can be directly taken from the experimental data in 
Ref.lH Using these data, we can completely determine 
the dynamics of any initial states using Eq. d3|. We 
then numerically find the optimal initial states pair, and 
the non-Markovianity of this system. The trace distance 
for the optimal initial states pair as a function of time 
obtained using our scheme is compared to the directly 
measured one in Fig. [I] As we see they are in excellent 
agreement. The non-Markovianity obtained from our 
method is 0.6, which are very close to the measured value 
0.58, both are in good agreement with the theoretical 
value 0.59. The error in our method is due to the error 
in the measurement of the 4 basis states. 


B. Open quantum walk system 

With the confidence of the method in two-dimension 
open system, we apply this method to a high-dimension 
open system which can not be reached previously. 
Here, we demonstrate our method to qualify the non- 
Markovianity of an one-dimension open quantum walk 
(QW) system. [22l f2~dl which is intrinsically a high- 
dimension system due to the ansatz coin. QW system is 
a generally interested system in quantum information, 
which has a lot of application in quantum computation. 
It has been shown that it is a nice tool to find new 
quantum algorithm and it can be used to constitute 
a universal model of quantum computation. [2^, [2(| 
In addition, QW has been experimentally realized in 
several different systems. [27H3ll ] The open QW has 
attracted a lot of attention recently. l22l 1211 

Here we study the discrete-time QW on a one- 
dimension lattice. The particle is located at one site at 
the beginning. At each step, it can move either to the left 
or to the right which is determined by the state of a coin: 
|L) (move left) or |i?) (move right). In quantum walk, the 
state of the coin can be a superposed state. Therefore, 
the state of the whole QW system (including the particle 
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FIG. 3: (Color online) The trace-distances as functions of 
time for the optimal state pairs obtained using our method 
for A'=0 (red line), 1 (blue line) and 2 (green line). 


in the lattice and the coin) is \i/j) = d ^x,d\x)\d) 
(in open QW, it should be a density matrix) where 
x = 0, ±1,±2, ••• are the location of the particle and 
d = L,R are the state of the coin. The operator to 
make up a single step of the QW can be defined as: 
W = TC , where T is the shift operator and defined as 

T = Ej\j - 1)01 ® |£><fl + 1 3 + 1)01 ® I R)(R\, 

C is the Hadamard coin operator defined as 
C = \L){L\ + \L)(R\ + \R)(L\ - \R)(R\). The 

state of the system after N step steps with the initial 
state |^o) can be obtained as W Nstep \ipo). 

The discrete-time QW can be implemented with single 
photon through an array of beam splitters in which the 
coin states are mimicked by the polarization degrees of 
freedom. For the open QW, the environment can be 
introduced through the coupling between the frequency 
and the polarization degrees of freedom of the photon 
similar to the method used in Ref. m In this case, the 
single step operator in the open QW can be modified as: 

= U 5t TC, (4) 

U 5t = [du e in ^|p)<p|®|w}<w|, (5) 

p—H,V 

where T and C are shift operator and Hadamard coin 
operator defined before. Ust couples the polarization)!/ 
or V) and environment to give non-Markovianity. n p 
is the index of refraction for different polarization state. 
5t p = is the time that operation TC takes place where 
L ~ 0.5 mm is the thickness of the beam splitters and v p 
is velocity of light with polarization p in the splitters. |w) 
is the environment state. For convenience, we take the 
environment as a delta function , that is |w) = — 

wi) + S(cu — W 2 )]. And uj\ = — wo, W 2 = fl + wp. 

We take the value of the parameters from the Ref. |l9|: 
wo = 7.2 x 10 12 s _1 , where A=780 nm, thus 

n = 2.4166 x 10 15 s -1 . n H = 1.554, n v = 1.545, St H = 
Lnji/c ~ 1.036 x 10~ n s and Sty ~ 1.030 x 10 -11 s. In 
addition, we use the periodic boundary condition for the 
system, i.e. \x m i n 1) = |x max ), |XmaxT 1) = |^mm) (see 


TABLE I: Comparing the number of initial states to measure 
between the present method and the direct scanning method 
to obtain the non-Markovianity of the quantum walk system. 
For the direct method, we assume 100 initial states to measure 
for each degree of freedom. The non-Markovianity of the 
system is calculated using the present method. 


number of initial states 

X Direc method Present method non-Markovianity 


0 

100 2 

4 

0.9512 

1 

100 6 

36 

0.9510 

2 

too 10 

100 

0.9428 


Fig.2), and we allow the system evolves for 20 steps. 

For convenience, we set the location of the particle 
x€[-X , X], then the total dimension of the system is 
N=2(2X + 1). To experimentally measure the non- 
Markovianity of this system, we need choose N 2 linearly 
independent initial states. As introduced before, we im¬ 
plement experiments to get the dynamics of the following 
initial states: 

1^) = k)M) (6) 

W = -^=(\xi)\d 1 )) + \x 2 )\d 2 )), (7) 

W = -^(\xi)\di)) + i\x 2 )\d 2 )) (8) 

here x,x\,X2=-X ■ ,X and d 1 d\ 1 d 2 =H, V. With the 
experimental data of these initial states, the dynamics 
of any initial state can be linearly constructed, and the 
optimal state pair can be found by numerically searching 
the configuration space. 

Because we do not have the experimental data at hand, 
we study the dynamics of the QW system numerically us¬ 
ing the exact diagonalization method. We first calculate 
the dynamics of the bases, then search the optimal initial 
states pair numerically described in previous paragraph, 
and determine the non-Markovianity. We collect the 
results of the trace distances of the optimal state pair for 
X=0, 1, 2 in Fig. [3] The non-Markovianity for X=0, 1, 
2 are given in Table Q] Actually, the result can be further 
confirmed via directly exact diagonalization method in 
which we obtain exactly the same results. 

We also compare the number of the initial states exper¬ 
imentally necessary to determine the non-Markovianity 
of the open QW system in Table Q] The number scales 
as 100 w , where N=4X+2 if previous direct scanning 
method is used. Clearly it is impossible to experimentally 
determine the non-Markovianity of the system for X >0 
using this method. However, by using our method, in 
which the number of initial states is only N 2 , we are able 
to measure the non-Markovianity of much larger systems. 
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IV. SUMMARY 

We have introduced a experimental method to quantify 
the non-Markovianity of high-dimension open quantum 
system. In our method, the scaling of the experiment is 
only N 2 which is dramatically reduced from exponential 
scaling of the conventional method. Therefore, the 
system which is intractable by the former method can 
be easily reached with the current method. After 
the benchmark with the well studied two-dimension 
open system, we demonstrate the method to the high 


dimension open quantum walk system. This method 
therefore opens up a new path to experimentally study 
the non-Markovianity of high dimension open quantum 
system, which was impossible previously. 
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